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Abstract 

Some new exact solutions of Einstein's field equations in a spatially homogeneous and anisotropic Bianchi 
type-V space-time with minimally interaction of perfect fluid and dark energy components have been ob- 
tained. To prevail the deterministic solution we choose the scale factor a{t) = Vi"e', which yields a time 
dependent deceleration parameter (DP), representing a model which generates a transition of the universe 
from the early decelerating phase to the recent accelerating phase. We find that for n > 1, the quintessence 
model is reproducible with present and expected future evolution of the universe. The other models (for 
n < 1), we observe the phantom scenario. The quintessence as well as phantom models approach to isotropy 
at late time. For different values of n, we can generate a class of physically viable DE models. The cosmic 
jerk parameter in our descended model is also found to be in good concordance with the recent data of 
astrophysical observations under appropriate condition. The physical and geometric properties of spatially 
homogeneous and anisotropic cosmological models are discussed. 
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1 Introduction 

Recent observations of type la supernovae suggest that the expansion of the universe is accelerating and two- 
thirds of the total energy density exists in a dark energy component with negative pressure [1] [2] (for a recent 
review, see Padmanabhan [3j; Copeland et al. [4]). In addition, measurements of the cosmic microwave back- 
ground [5] and the galaxy power spectrum [B] also indicate the existence of the dark energy. However, the 
observational data are far from being complete (for a recent review, see Perivolaropoulos [7]; Jassal et al.[S]). 
It is not even known what is the current value of the dark energy effective equation of state (EoS) parameter 
which lies close to —1: it could be equal to —1 (standard ACDM cosmology), a little bit 
upper than —1 (the quintessence dark energy) or less than —1 (phantom dark energy). While the possibility 
a; ^ — 1 is ruled out by current cosmological data from SN la (Supernovae Legacy Survey, Gold sample of 
Hubble Space Telescope) [TO], CMB (WMAP, BOOMERANGE) [H [12] and large scale structure (Sloan 
Digital Sky Survey) dataFT5|, the dynamically evolving DE crossing the phantom divide line (PDL) (w = —1) is 
mildly favoured. The simplest candidate for the dark energy is a cosmological constant A, which has pressure 
Specifically, a reliable model should explain why the present amount of the dark energy is so 
small compared with the fundamental scale (fine-tuning problem) and why it is comparable with the critical 
density today (coincidence problem) [?|. That is why, the different forms of dynamically changing DE with 
an effective equation of state (EoS), 

I p{d-'^) < —1/3, have been proposed in the literature. Some 
other limits obtained from observational results coming from SNe la data [13] and combination of SNe la data 
with CMBR anisotropy and galaxy clustering statistics |B] are —1.67 < w < —0.62 and —1.33 < w < —0.79, 
respectively. The latest results in 2009, obtained after a combination of cosmological datasets coming from CMB 
anisotropics, luminosity distances of high redshift type la supernovae and galaxy clustering, constrain the dark 
energy EoS to -1.44 < w < -0.92 at 68% confidence level 



Moreover, in recent years Bianchi universes have been gaining an increasing interest of observational cos- 
mology, since the WMAP data [IS] — [13 seem to require an addition to the standard cosmological model with 
positive cosmological constant that resembles the Bianchi morphology [H] — [23]. According to this, the universe 
should achieved a slightly anisotropic special geometry in spite of the inflation, contrary to generic inflationary 
models and that might be indicating a non-trivial isotropization history of universe due to the presence of an 
anisotropic energy source. In principle, once the metric is generalized to Bianchi types, the EoS parameter of 
the fluid can also be generalized in a way conveniently to wield anisotropy with the considered metric. In such 
models, where both the metric and EoS parameter of the fluid are allowed to exhibit an anisotropic character, 
the universe can exhibit non-trivial isotropization histories and it can be examined whether the metric and/or 
the EoS parameter of fluid evolve toward isotropy. Thus, the Bianchi models which remains anisotropic are 
of rather academical interest. The study of Bianchi type V cosmological models create more interest as these 
models contain isotropic special cases and permit arbitrary small anisotropy levels at some instant of cosmic 
time. This property makes them suitable as model of our universe. The homogeneous and isotropic Friedman- 
Robertson- Walker (FRW) cosmological models, which are used to describe standard cosmological models, are 
particular case of Bianchi type /, V and IX universes, according to whether the constant curvature of the 
physical three-space, t — constant, is zero, negative or positive. 

Recently, Yadav [H] and Kumar & Yadav [2S] have dealed with a spatially homogeneous and anisotropic 
Bianchi type-V DE models by considering constant DP. In this paper, we have considered minimally interacting 
perfect fluid and DE components with time dependent DP within the framework of totally anisotropic Bianchi-V 
models. This paper is organized as follows: the metric and the field equations are presented in Sect. 2. Sect. 
3 deals with the exact solutions of the field equations and the physical behaviour of the model. Section 4 deals 
with cosmic jerk parameter. Discussions and concluding remarks are given in Section 5. 

2 The Metric and the Field Equations 

We consider the space time metric of the spatially homogeneous and anisotropic Bianchi type-V of the form 

ds^ ^ ^dt^ + A^dx^ + e^""' [B^dy^ + C^dz^] , (1) 
where the metric potentials A, B and C are functions of cosmic time t alone and a is a constant. 

We define the following physical and geometric parameters to be used in formulating the law and further in 
solving the Einstein's field equations for the metric (jlj. 

The average scale factor a of Bianchi typc-V model ([T]) is defined as 

a=(^BC)i (2) 

A volume scale factor V is given by 

V = a^ = ABC. (3) 
We define the generalized mean Bubble's parameter H as 

I 
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H ^ -{H,+Hy + H,), (4) 



where Hx = ^, Hy = ^ and = ^ are the directional Bubble's parameters in the directions of x, y and z 
respectively. A dot stands for differentiation with respect to cosmic time t. 

From Eqs. ©-(HI), we obtain 

IV h I ( A B C\ 
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The physical quantities of observational interest in cosmology i.e. the expansion scalar 0, the average anisotropy 
parameter Am and the shear scalar ct^ are defined as 
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where AiJ^ = Hi — H{i = x, y, z) represents the directional Hubble parameter in the direction of x, y, z respec- 
tively. Am = corresponds to isotropic expansion. 



The Einstein's field equations ( in gravitational units SttG = c = 1) read as 



1 



(de)i 



(9) 



where T^J™-** and yj'^'^-'* are the energy momentum tensors of perfect fluid and DE, respectively. These are given 

by 



= diaghp('"\p('"),p('"),p(™)], 
= diaghl,w(™\cj('"\w(™']p'", 



(10) 



and 



rf^)'=diag[-p('^'^),p('''=),p('^^),p('^'^)], 



= diag[-l,a;^'"=^cJ^''^^cJ^'"=>^'"=^ (11) 

where p^™^ and are, respectively the energy density and pressure of the perfect fluid component or ordinary 
baryonic matter while / p{m) is its EoS parameter. Similarly, p^'^'^^ and p'^'^'^'^ are, respectively the 

energy density and pressure of the DE component while I p^"*"^) is the corresponding EoS parameter. 

We assume the four velocity vector = (1,0,0,0) satisfying v}u.j — —1. 

In a co-moving coordinate system (u' = Jg), Einstein's field equations (|9|) with (fTO|) and (fTT|) for B-V metric 
([J) subsequently lead to the following system of equations: 
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The law of energy-conservation equation (T.j' = 0) yields 

p{m) ^ ^ a;(™))p(™)i7 + p^'^'^^ 3(1 uj^^'''^)p^'^'''^ H = 0. 
The Raychaudhuri equation is found to be 

^? = -(l + 3c^(''^))p(''^)-i02_2a2. 



(12) 

(13) 
(14) 
(15) 
(16) 
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3 Solution of the Field Equations and its Physical Significance 

In order to solve the field equations completely, we first assume that the perfect fluid and DE components 
interact minimally. Therefore, the energy momentum tensors of the two sources may be conserved separately. 

The energy-conservation equation (y^™^*^ — Q) of the perfect fluid gives 

p('")+3p('")(l+cj('"))iJ = 0, (19) 

where as energy-conservation equation (X'^''"^)*-' — Q) of the DE component leads to 



P 



(de) 



+ 3p('*^)(a;('*") + = 0. (20) 



Following, Akarsu and Kilinc (2010), Yadav (2011) and Kumar & Yadav (2011), we assume that the EoS 
parameter of the perfect fluid to be constant, that is, 

"^"^ = 7^ - (21) 

while w*^'^'^-' has been permitted to be a fimction of time since the current cosmological data from SN la, CMB 
and large scale structures mildly favour dynamically evolving DE crossing the PDL as discussed in previous 
section. 



Eq. (|T9|) can be integrated to obtain 

= poa"^^""*'^^ (22) 

where po is a positive constant of integration. 



Integrating (|16p and engrossing the constant of integration in B or C, without any loss of generality, we obtain 

= BC. (23) 



We now use the technique following Kumar and Yadav 25: to solve Einstein's field equations ([T2l) — (fT5|) . 
Subtracting ([T^ from (IT^ , (IT^ from ([T^ , and from P3|) and taking second integral of each, we obtain the 
following three relations respectively: 

A . f . f df 
B 



di exp /ci / — , (24) 



^ = d,e^p(k,J (25) 

and 

^ = da exp (^/C3 J , (26) 

where di, d2, ds, ki, k2 and fcs are constants of integration. From ([M)) — ([26 |) and (|23p . the metric functions can 
be explicitly obtained as 

Ait) = a, (27) 

B{t) = ma exp ( £ / ^ ) , (28) 



a 



dt 



C{t) = -e^p[-i -], (29) 

W^), £=i^i±M, d2^d-^\ k2 = -k,. (30) 

Finally, following Saha et al. [26] . we take following ansatz for the scale factor, where increase in term of time 

evolution is 

a{t) = (31) 



where 

m 
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where n is a positive constant. This ansatz generalized the one proposed by Amirhashchi et al. [37]. In 
hterature it is common to use a constant deceleration parameter [231 Hi] , [IB] — [32] as it duly gives a power law 
for metric function or corresponding quantity. The motivation to choose such time dependent DP is behind 
the fact that the universe is accelerated expansion at present as observed in recent observations of Type la 
supernova [T] [21 [TOl 1231 133] and CMB anisotropies [3S] — [3Z] and decelerated expansion in the past. Also, the 
transition redshift from deceleration expansion to accelerated expansion is about 0.5. Now for a Universe which 
was decelerating in past and accelerating at the present time, the DP must show signature flipping [55] — [ID]. 
So, there is no scope for a constant DP at the present epoch. So, in general, the DP is not a constant but time 
variable. The motivation to choose such scale factor (|3T]) yields a time dependent DP. 
Using pil) into (|?7)) - (|3n)) . we get the following expressions for scale factors: 




t 



n=Q.5 — - — n = \ n=2 



Figure 1: The anisotropic parameter Am versus t. Here i — 1 



A(i) = (re*)5, (32) 
B{t) = m(i"e*)^ exp f {re'yidt] , (33) 



C{t) = m-'{re')^ exp l^-i J (ee'y^dtj. (34) 

The expressions for physical parameters such as directional Hubble parameters {H^, Hy, Hz), the Hubble 
parameter {H), scalar of expansion (0), shear scalar (ct), spatial volume V and the anisotropy parameter {Am) 
are, respectively, given by 
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Figure 2: The deceleration parameter q versus t 



a^^f{t^ey\ (39) 
F= (re*)^exp(2ax), (40) 

A. = f + (41) 

It is observed that at i = 0, the spatial volume vanishes and other parameters 9, a, H diverge. Hence the model 
starts with a big bang singularity at < = 0. This is a Point Type singularity (MacCallum 1971) since directional 
scale factor A{t), B{t) and C{t) vanish at initial time. 

Figure 1 depicts the variation of anisotropic parameter {Am) versus cosmic time t. From the figure, we ob- 
serve that Am decreases with time and tends to zero as f — > c» for all values of n. Thus, the observed isotropy 
of the universe can be achieved in our derived model at present epoch. The shear tensor also tends to zero in 
this model. 



We define the deceleration parameter q as usual, i.e. 



Using Eq. (|3Tj) into Eq. (|42j), we find 



aa a , , 

1 = -— = —FH- 42 



[n + ty 



From Eq. (j43|) . we observe that g > for t < v2ri — n and g < for < > \/2n — n. It is observed that for 
< n < 2, our model is evolving from deceleration phase to acceleration phase. Also, recent observations of 
SNc la, expose that the present universe is accelerating and the value of DP lies to some place in the range 
— 1 < g < 0. It follows that in our derived model, one can choose the value of DP consistent with the observation. 
Figure 2 graphs the deceleration parameter (q) versus time which gives the behaviour of q from decelerating to 
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Figure 3: The plot of DE EoS parameter oj'^'^") versus t. Here po = 1, a = 1, ^ = 1, = 0.5 



accelerating phase for different values of n. 

The energy density (p*-™') of perfect fluid, the pressure 

(p('ie)) Qf component, DE density {p'^'^")) and EoS 

parameter (a;^'*'^^) of DE, for this model are given by 

=po(t"e*)"i(i+"'"'), (44) 
pi''-) = _^ + l)' + _ ^(re*)-3 + a2(<»e*)-i - c^(™)po(t"e*)-i(i+"""'), (45) 
pi''-) = -po(re*)~5(i+'^""') + ^ ^ l) V £(re*)-3 - 3a\ee')-\ (46) 
^(,,) ^ I (f + 1)' - F + - «^(t"e*)-i + c.("')po(t"e*)-i(i+"""') 

Figure 3 depicts the variation of DE EoS parameter uji''-) versus cosmic time t. We observe from the figure that 
for n < 1, uji''-) varies from non-dark region crossing the PDL (uji''-) — —1) and ultimately approaches to the 
phantom region [uji''-) < —1). But for n > 1, the variation of uji''-) starts from cosmological constant region 
(ujid-) = -1) and finally approached to the quintessence region (o;'^'*'^^ > — 1). Therefore, we observe that for 
n > 1, the variation of w*^'*'^-' in our derived model is consistent with the recent observations of SNe la data [14, . 
SNe la data with CMBR anisotropy and galaxy clustering statistics [B]. 

The dark energy with {jji''-) < —1, the phantom component of the universe, leads to uncommon cosmologi- 
cal scenarios as it was pointed out by Caldwell et al. [42]. First of all, there is a violation of the dominant 
energy condition (DEC), since p + p < 0. The energy density grows up to infinity in a finite time, which leads 
to a big rip, characterized by a scale factor blowing up in this finite time. These sudden future singularities 
are, nevertheless, not necessarily produced by a fluid violating DEC. Cosmological solutions for phantom matter 
which violates the weak energy condition were found by Dabrowski et al. [43]. Caldwell [44], Srivastava [45] , 
Yadav [21] have investigated phantom models with uji'''^) < — 1 and also suggested that at late time, phantom 
energy has appeared as a potential DE candidate which violets the weak as well as strong energy condition. 

The left hand side of energy conditions have been depicted in Figures 4 and 5 for different values of n. From 
Figure 4, for n = 0.5 (i.e phantom model) (see. Figure 3) , we observe that 

(i) pi"^") > 0, (m) p^'^^^ + pi"^"^ < 0, {tit) p'^'^') + ip'^^"^ < 0. 
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Figure 4: The plot of energy conditions versus i for n = 0.5. Here po — I, i — a — 1, w^™^ — 0.5 



Thus, from above expressions, we observe the phantom model violates both the strong and weak energy condi- 
tions, as expected. 

Furthe, from Figure 5, for n > 1 (i.e quintessence model) (see, Figure 3), we observe that 
(z) P^'^"^ > 0, (m) p'-'^"'^ > 0, (zm) p^'^"'^ + Sp^"^"^ < 0. 

Thus, the quintessence model violates the strong energy condition as the same is predicted by current astro- 
nomical observations. 



The perfect fluid density parameter (ri*^™^) and DE density parameter (51^'^'^^) are given by 

f]M = 4^o(y + l)'(tV)-i(i+"<'"'). 



(48) 
(49) 



Thus the over all density parameter (fi) is obtained as 

Figure 6 depicts the variation of the density parameter (fl) versus cosmic time t for different values of n during 
the evolution of the universe. From the Figure 6, it can be seen that the total energy density ft tends to 1 for 
sufficiently large time which is supported by the current observations. 



(50) 



4 Cosmic Jerk Parameter 

A convenient method to describe models close to A CDM is based on the cosmic jerk parameter j, a dimensionless 
third derivative of the scale factor with respect to the cosmic time [36] — [50]. A deceleration-to-acceleration 
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Figure 5: The plot of energy conditions versus t for n — \. Here pq = 1, I = a — \, w^™^ = 0.5 



transition occurs for models with a positive value of jo and negative go- Flat A CDM models have a constant 
jerk J = 1. The jerk parameter in cosmology is defined as the dimensionless third derivative of the scale factor 
with respect to cosmic time 

1 a 



and in terms of the scale factor to cosmic time 



ja^H^)" 
2H^ 



(51) 



(52) 



where the 'dots' and 'primes' denote derivatives with respect to cosmic time and scale factor, respectively. One 
can rewrite Eq. ([5T|) as 



H 



Eqs. (|43|) and dSS]) reduce to 



Jit) = 1 



6n 



(53) 



(54) 



{n + tf (n + l)3- 

This value overlaps with the value j ~ 2.16 obtained from the combination of three kinematical data sets: the 
gold sample of type la supernovae [5], the SNIa data from the SNLS project [TD], and the X-ray galaxy cluster 
distance measurements [51] for 

(55) 



t = 3.45x lO-^A-^-n, 



where A = 104n[8.41 + 1.45 \/(14.4n + 33.6)] i 



5 Discussion and Concluding Remarks 

In this paper, we have studied a spatially homogeneous and anisotropic Bianchi type-V space time filled with 
perfect fluid and anisotropic DE possessing dynamic energy density. The field equations have been solved exactly 
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Figure 6: The total energy density parameter D. versus t. Here a — I — 1 



with suitable physical assumptions. The solutions satisfy the energy conservation Eq. (|17|) and the Raychaud- 
huri Eq. (fT8|) identically. Therefore, exact and physically viable Bianchi type-V model has been obtained. It is 
to be noted that our procedure of solving the field equations is altogether different from what Kumar & Yadav 
[25] have adopted. Kumar & Yadav [35] have solved the field equations by considering the constant DP whereas 
we have considered time dependent DP. As we have already discussed in previous Sect. 3 that for a Universe 
which was decelerating in past and accelerating at the present time, the DP must show signature flipping (see 
the Refs. Padmanabhan and Roychowdhury [35], Amendola [321, Riess et al. [10] and so, there is no scope for 
a constant DP. The main features of the model are as follows: 

• For different values of n the anisotropic parameter Am tends to zero for sufficiently large time. Hence, the 
present model is isotropic at late time which is consistent to the current observations. 

• The present DE model has a transition of the universe from the early deceleration phase to the recent 
acceleration phase (see, Figure 2) which is in good agreement with recent observations [52|. 

• In the present study we find that for n > 1, the quintessence model is consistent with present and expected 
future evolution of the universe. The quintessence model approaches to isotropy at late time (see. Figure 1 & 
3). The other models (for n < 1), we observe the phantom scenario. 

• The derived phantom model violates both the strong and weak energy conditions whereas the quintessence 
model violates only the strong energy condition (see. Figures 4 & 5). 

• The total density parameter (f2) approaches to 1 for sufficiently large time (see. Figure 6) which is repro- 
ducible with current observations. 

• The cosmic jerk parameter in our descended model is also found to be in good agreement with the recent 
data of astrophysical observations namely the gold sample of type la supernovae [9], the SNIa data from the 
SNLS project [10], and the X-ray galaxy cluster distance measurements [STj . 
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• Our special choice of scale factor yields a time dependent deceleration parameter which represents a model 
of the Universe which evolves from decelerating phase to an accelerating phase whereas in Yadav [24) . Kumar 
& Yadav [25^ only the evolution takes place either in an accelerating or a decelerating phase. 

• For different choice of n, we can generate a class of DE models in Bianchi type-V space-time. It is observed 
that such DE models are also in good harmony with current observations. Thus, the solutions demonstrated in 
this paper may be useful for better understanding of the characteristic of anisotropic DE in the evolution of the 
universe within the framework of Bianchi type-V space-time. 
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